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A theoretical scheme for an experimental implementation involving bisolitonic matter 
waves from an attractive Bose-Einstein condensate, is considered within the framework 
of a non-perturbative approach to the associate Gross-Pitaevskii equation. The model 
^ ' consists of a single condensate subjected to an expulsive harmonic potential creating 

^\ ' a double-condensate structure, and a gravitational potential that induces atomic ex- 

^SJ , changes between the two overlapping post condensates. Using a non-isospectral scatter- 

ff^ ■ ing transform method, exact expressions for the bright-matter-wave bisolitons are found 

^^ ' in terms of double-lump envelopes with the co-propagating pulses displaying more or 

1^* ' less pronounced differences in their widths and tails depending on the mass of atoms 

^ - , composing the condensate. 

s«.i^ ■ Keywords: Bose-Einstein condensate; gravity; anti trap; Gross-Pitaevskii equation; non- 

isospectral scattering transform; matter-wave bisolitons. 

>: 

j^ ' 1. Introduction 

H 

Cu ' Outstanding progress has been made in the domain of laser physics thanks to the 

discovery of Bose-Einstein condensation L^J^, a phenomenon observed in Bose liq- 
uids upon cooling and preparation in optical confinements'^. Atom lasers' 'I 'thus 
consist of coherent matter-wave structures that recently turne d to be ideal inputs 
in quantum interferometry metrology such as the double-slit i^^ and Michelson- 
Morley 1^ interferometries. Given the robustness of their solitonic shapes, matter- 
wave lasers also represent an enormous potential for holographic applications involv- 
ing coherent beams of optical fields which are non-dest ructiv e under interferences, 
wave mixings and many-body recombination processes '^ ' I 

The present work addresses a particularly ineresting matter-wave laser implemen- 
tation in which the input atom laser is a double-lump pulse, obtained by splitting 
a single Bose-Einstein condensate(BEC) into two parts. Such implementation has 
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relevant application in the double-slit holography '^^ where two overlapping con- 
densates are created by splitting one single condensate, e.g. by a far-ofF-resonance 
laser bar rier l^^, leading to two atomic populations in nearly degenerate hyperfine 
levels LUJ. In many past theoretical studies devoted to the last subject, the far-ofF- 
resonance splitting is described by a double- well potential ^^. A main motivation 
for the choice of a double- well potential is the presence of a potential hump that cuts 
the condensate into two slices, and of two minima allowing harmonic confinement 
of each of the two overlapping condensates. However, in a strict sense the mean- 
field picture predicts that the governing mechanism of formation of matter-wave 
solitons is the interparticle interaction and the trap interefers as a post process, 
the role of which is to confine matter- wave solitons thus permitting full control and 
experimental manipulations of the resulting laser beam. Therefore, in studying the 
far-off-resonance splitting process of a single BEC it is enough considering only the 
hump part of the double-well potential, a pa rt which can well be represented by an 
harmonic expulsive (i.e. anti-trap) potential^ I I I I 

Matter-wave lasers are relatively more sensitive to a gravitational field than the 
classic photon laser for the matter-wave laser source is an excited population of 
atoms. Particularly in the case of condensates consisting of vapor atoms the free fall 
stands for a relevant process promoting atomic exchanges between two overlapping 
condensates, and this process is expected to be dominant over scattering processes 
such as thermal feedings and many-body recombinations in the case of heavy atoms 
BECs. From the theoretical standpoint, the influence of gravity o n Bos e systems 
has been discussed, flrstly in the context of non-interacting particles ^ I I and more 
recently in the case of a system of weakly interacting Bose particles ^^. In this last 
work, following the WKB approximation the authors obtained that the condensate 
wavefunction could be approximated by solutions of a perturbed linear Schrodinger 
equation. On the basis of this assumption, they derived expressions of the conden- 
sate wavefunction as combinations of Bessel functions. 

In this work we shall develop a non-perturbative theory for the mean-field dynamics 
of an atomic population consisting of two condensates held together by an harmonic 
anti-trap potential and exchanging atoms by means of the gravitational field. As 
already indicated, a BEC phenomenon of this kind has established physical impor- 
tance l^^^-^^!-^^^^^^^^^ and it will emerge below that in the specific context of BEC 
systems with relatively heavy atoms, the gravity can strongly affect both shape 
profiles and th e dynamics of pulses in the matter-wave bisoliton created by the 
anti-trap l^^l^ so that perturbative techniques, including variational methods, are 
very limited in their description. 



2. Model and analytical bright-matter bisoliton solutions 

Consider the following Gross-Pitaevskii equation, which governs the groundstate 
wavefunction of a BEC system in an anti-trap potential and subjected to the grav- 
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P,z-m^z2_2^|^|2 ^^ (^) 



where ^ is the (positive) mean-field interatomic interaction coefficient related to the 
s — wave scattering length, m is the atomic mass and Pa = rng the associate weight 
with g the acceleration of gravity. u}\\ is the characteristic frequency for atomic os- 
cillations along the axis of graviational force, as we are interested in the dynamics 
of the above system along the last axis l^ we shall assume V^ = d^/dz^. 
To gain the best formulation of the dynamics of BEC systems described by the 
mean- field Eq. ([T]), we follow the Inverse-Scattering- Transform (1ST) method which 
offers an interesting non-perturbative framework for the study of exactly integrable 
nonlinear partial differential equations. To this last end, remark to start that when 
the gravity and the expulsive potential are absent Eq. ([T]) reduces to the homo- 
geneous Nonlinear Schrodinger (NLS) equation which, within the framework of 
the 1ST method, gives rise to an isospectral Lax-pair-type eigenvalue problem with 
universal conservation laws . However, in the presence of perturbation, this equa- 
tion can become a rather complicated mathematical problem and is rarely exactly 
integrable. Nevertheless, exploiting rece nt developnients i n the treatment of non- 
isospectral Lax-pair eigenvalue problems l ^'MoMtJ | oU | ol | -J^ I^ non- universal conserva- 
tion laws can be derived for a certain class of perturbations such that the associate 
NLS equations are exactly integrable. In the specific case of Eq. ffl such conser- 
vation laws have been constructed in the cases Pa y^ but aj|| = , Pa ~ but 
wii 7^ 1^, Pa ^ and a;|| ^ Lt^. Following the method for this last case, the two 
pairs of linear eigenvalue problems associate with the perturbed NLS Eq. Q are 
defined as: 

U2 + i Xu ^ ip{z, t) V, 

Vz - iXv ^ -tl;*{z,t)u, (2) 

Ut ~ Au — B V = 0, 

vt-Cu + Av = 0, (3) 

where A , B, C are coefficients of the matrix associate with the second linear sys- 
tem dl]) ''^^ 1 ^ I, u{z,t) and v{z,t) are space-time dependent functions forming el- 
ements of a two-component wavevector while A is the associate eigenvalue. The 
quantity ijj{x,t) is the solution wanted , but now plays the role of IST's scattering 
potential. According to the general rule I I I I of the 1ST technique, the appropri- 
ate space-time evolutions of the matrix elements A, B and C should be determined 
unambiguously provided conservation laws for the set ^ and ([3]) are explicitely 
formulated. These conservation laws follow from the condition of compatibility be- 
tween ^ and dS]) corresponding to a derivation of © with respect to t and ^ 
with respect to z and comparing the four resulting equations. This leads to the 
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following set of three coupled evolution equations: 

A,+tXt = -Cij{z,t)-B4,*{z,t),, (4) 

(In B), +2iX = il)t{z, t)/B - 2A ^{z, t)/B, (5) 

(InC)^ -2i\ = -i^tiz, t)/C - 2A il;*{z, t)IC, (6) 

where the eigenvalue A is assumed space-time dependent as a result of the pertur- 
bation in the NLS (H]), hence the non-isospectral character of the 1ST eigenvalue 
problem. 

As we are interested in the matter-wave structures with pulse shape, we need to 
operate a judicious choice of the scattering potential. For the shape profile of cur- 
rent interest the scattering potential must vanish asymptotically as z — > ±cx) at any 
propagation time t, meaning that a reflectionless potential stands for an appropri- 
ate choice. In the specific case of the perturbed NLS Eq. ([1]), a best representation 
for the initial soliton solution with pulse shape is: 



V'(z, i = 0) = -2A(z, t = 0) sech 



\{z',t^O)dz' 



(7) 



where X{z,t = 0), the eigenvalue of the 1ST problem at the initial time t = 0, is 
assumed to vary in both space and time in accordance with the non-isospectral 
character of the associate eigenvalue problem. To put the last remark on a firm 
ground, it suffices to notice that by combining the compatibility relations (H])-® 
with the generic equation ((T)) , we can derive a single compatibility equation in terms 
of the sole eigenvalue A i.e.: 

V^ - Xt ~ 2{iX)l = 0, V{z)^ Paz-mup"^. (8) 

This compatibility equation indeed gives clear evidence of a dependence of A on 
the space and time coordinates, establishing thus the non-universal character of 
spectral parameters of the eigenvalue problems ©-([S]). It turns out that A can be 
written explicitely as A(z,i), and in general any solution provided by ([8]) must be 
generic of eigenvalue(s) of the 1ST scattering potential defined in ([7]). 
Now groupi ng all the preceeding key considerations in the standard steps of the 
1ST method U4|zd|.j4| ^ ^^ arrive at the following exact one-soliton solution for the 
perturbed NLS Eq. (P): 

'ipiz,t) ~ ipQ{z,t)sech[2^{z,t)] exp[— 193(2, t)], 

where 



M^^,t) = ^oG{z)F,{t), (9) 



^z,t)^-^F.it) 



(fiiz^t) 



2^/2m^ 




^-^ f^^[{F\t')]dt\ 



^[{F^{t')]dt', 



(10) 



(11) 



(12) 
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F{t) = Fi{t) + iFr{t) with 

^1 (l + A^)tanh(Kt) ^ ^ sech^JKt) 
"^ ' X l + XHixnh^{Kt)' '^ ' l + XHixnh'iKt)' ^ ' 

and 

K = -2iLj||v/m//i, A = 1/k. (14) 

Note that the quantity -d^ in © depends entirely on amplitudes of the conjugate 
pair of eigenfunctions (i.e. Jost functions) of the 1ST scattering problem L!^, hence 
it does not depends on space and time and can therefore readily be regarded as an 
arbitrary amplitude of the 1ST one-soliton solution. 

3. Numerical analysis of the analytical bisoliton solution 

Consider an hypothetical Bose- Einstein condensation experiment aimed at creating 
an atom laser for a double-slit interferometer. Starting from an initial setup consist- 
ing of one single BEC with identical masses and an attractive mean- field interaction 
between the atoms, we switch on adiabatically an harmonic anti-trap potential to 
produce a smooth barrier hump (i.e. a non-flat potential barrier) that can cut the 
single BEC into two slices. This gives rise to two overlapping BECs, and due to 
the presence of a finite attractive mean-field interaction in the two condensates we 
generate twinned matter-wave pulses of equal amplitudes. Next, suppose that the 
created bisolitonic matter-wave structure is subjected to the gravity force, so that 
the two constituent pulses are accelerated one with respect to the other. Overall, 
the gravity potential will cause breaking of symmetry of the harmonic expulsive 
potential and from the viewpoint of the Gross-Pitaevskii Eq. ([T|) , we expect a di- 
rect incidence of this symmetry breaking on the matter-waves structures. Our goal 
in this section is to examine what the analytical solution to the Gross-Pitaevskii 
equation for the system just described, obtained in ([9]) via the 1ST technique, pre- 
dicts in terms of the influence of combined effects of the expulsive potential and 
gravity on the dynamics of the double condensate at different representative orders 
of magnitude of the atomic mass m. 

First, we start with an analysis of the effects of the gravitational force on spatial 
profiles of the matter-wave bisoliton. Fig. ([T]) displays the magnitude 1-0(^,0)1^ of 
the matter- wave intensity as given by ([9]) at the initial time i = 0, for a fixed value 
of the amplitude vnuil of the harmonic barrier (i.e. mwn = 0.5) but different atomic 
weights. More precisely, from the first to the last graph values of P^ (in arbitrary 
units) are 0, 0.01, 0.1, 0.4, 0.5 and 0.8 respectively. The first graph corresponds 
to a context where the gravity force is absent (i.e. Pa = 0) and only the anti-trap 
potential is switched on. As one sees, the two pulses in the bisoliton structure ^ 
form a pair of perfect twins, characterized by equal tails and widths and hence 
should carry same amount of energy. Now when we also switch on the gravitational 
field, the perfect twinning of pulses breaks down as it emerges through shapes of 
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Fig. 1. Amplitude |i/i(2,0)| of tiic initial matter-wave intensity versus position z. For all the six 
graphs muj?^ = 0.5 and Pa is varied as foUowrs, from the top left to the last graph: 0, 0.01, 0.1, 0.4, 
0.5, 0.8. 



the matter- wave intensity in the five last graphs of fig. ([ij . Particularly remarkable 
is the difference of amplitudes among the two pulses, which is more and more pro- 
nounced as the atomic weight increases. At sufficiently large values of Pa, one of 
the twin pulses seems to collapse reflecting a strong infiunce of the gravity tend- 
ing to confine all atoms within one of the two condensates. The last behaviour is 
sam ewh at comparable with the well-known phenomenon of macroscopic self trap- 
ping '^^, but except that in the present context atoms are exchanged between the 
two overlapping condensates by means of classical motion while in the common 
macroscopic self-trapping phenomenon, the density of population imbalanceis de- 
termined by quantum tunnelings via a Josephson-like phase dynamics I ^ ^ [ It 
should equally be stressed that in real experiments, the account of gravity must 
favour non-symmetric positions of the two overlapping condensates with respect to 
a common axis coinciding with the axis of the gravity force. In such context the free 
fall process provides a relevant enhancing mehanism for atomic exchanges between 
the two condensates for heavy atoms that are expelled from an upper condensate 
to a down condensate by means of thermal ejections (i.e. thermal evaporation). 
We also examined profiles of the moving twinned pair, still for different values of 
the atomic weight and keeping mwn = 0.5. Since the integrals with respect to time 
in (fTO| and (|11|) are not tractabe analytically, we applied a numerical integration 
scheme assuming i = as the initial time. Figs. ^ and ^ summarize variations of 
the matter-wave intensity in space at different times, for different atomic weights. 
On fig. ([2]), the left column displays spatial profiles of the bisoliton intensity for 
Pa = at five different times namely t = 0, 1.5, 2.5, 3.5 and 4 from the top graph 
respectively, while the right column are spatial profiles at same times but for an 
atomic weight Pa = 0.01. Fig. ([3]) is a continuation of fig. 1^ where the left col- 
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Fig. 2. Spatial profile of the intensity \ip{z,t)\ of the matter- wave bisoliton for Pa = (left column) 
and Pa = 0.01 (right column). From the top to the bottom graphs times are choosen as follows, 
in arbitrary units: t = 0, 1.5, 2.5, 3.5, 4 
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Fig. 3. Spatial profile of the intensity |^(2,t)| of the matter-wave bisoliton for Pa = 0.05 (left 
column) and Pa = 0.1 (right column). From the top to the bottom graphs times are choosen as 
follows, in arbitrary units: t = 0, 1.5, 2.5, 3.5, 4 
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umn corresponds to Pa = 0.05 and the right cohimn to Pa — 0.1. According to 
curves of the two figures, in general the twin pulses composing the matter-wave 
bisoliton described by the analytical solution ([9]) propagate in opposite directions. 
As the graphs for finite weights indicate, when Pa increases the twin pulses aquire 
different amplitudes while still moving in opposite directions. At sufficiently long 
propagation time, their amplitudes both decrease and the bisoliton tends to recover 
its perfect twin structure. 

An exact explanation of the apparent recovery of perfect twinning of the two pulses 
in the bisoliton after a long propagation time is not evident, as this behaviour 
emerges through numerical simulations of the full expression (j9|) at finite times and 
positions. However, we suspect this behaviour to originate from a complex variation 
of the prefactor ^o{z,t) when both time and space are simultaneously nonzero. In 
fact, the matter- wave amplitude is actually not constant but changes with space 
and time suggesting a non-topological feature for the bisoliton solution found in (j9|). 

4. Conclusion 

In summary we developed a non-perturbative approach to the Gross-Pitacvskii 
equation in the presence of an harmonic non-symmetric expulsive potential, in or- 
der to gain a more consistent insight about the dynamics of BEG systems under 

oo no] 

combined effects of an expulsive potential and the gravitational field ^^'^ We 
found that the expulsive potential, which has shape of a smooth space- varying (i.e. 
harmonic) barrier, favors twinned matter-wave structures and thus acts like an 
atomic beam splitter. The gravity was found to break the symmetry of this har- 
monic barrier turning the perfect twin-pulse bisoliton into an asymmetric bisoliton. 
This asymmetry emerged through their distinct tails, widths and positions with 
respect to a common initial position along the gravity axis, and appeared to be 
more and more pronounced as the atomic weight was increased. However, while the 
spatial profile of bisoliton suggests a monotonic change in initial amplitudes of the 
two pulses, the space-time evolution instead reveals a complex behaviour involving 
possible recovery of the perfect pulse twinning after a sufficiently long propagation 
time. This behaviour has been infered to a non-topological (i.e. dynamical) feature 
of the bisoliton. 

Dynamical properties of BECs under combination of the expulsive harmonic barrier 
and the gravity that emerge in the present study, reveal a quite rich physics some 
aspects of which have already be experimentally predicted I ^^JMl | ^^I^^M4 | /io | r^^^ 
mathematical method adopted, namely the 1ST technique, clearly provides a frame- 
work where this rich physics can be captured far better than within perturbative 
approaches, or some variational treatments of the anti-trap and/or gravity. 
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